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We review the renormalization of the electroweak sector of the standard model. 
The derivation also appUes to the minimal supersymmetric standard model. We 
discuss regularization, and the relation between the threshold corrections and the 
renormalization group equations. We consider the corrections to many precision 
observables, including Mw and siv? 0^^ . We show that global fits to the data 
exclude regions of supersymmetric model parameter space and lead to lower bounds 
on superpartner masses. 



1 Introduction 

1.1 The standard model and beyond 

Over the past decade we have seen the standard model (SM) tested with in- 
creasing precision by a large number of experiments. It is quite remarkable 
that as the measurements become more accurate and encompass more observ- 
ables, the standard model is repeatedly confirmed to greater precision. When 
discrepancies are found, they always seem to evaporate with subsequent data 
and/or analyses. From the point of view of precision measurements there is no 
need for physics beyond the standard model. 

With so much direct experimental evidence confirming the standard model 
(and no evidence of discrepanciesu) it might appear somewhat mysterious that 
so much time and effort is spent studying models beyond the standard model. 
There are many reasons why people believe the standard model with a light 
Higgs boson is not the whole story. Many agree that new physics near the 
electroweak scale is necessary in order to resolve the hierarchy problem. They 
point out that in order for the huge hierarchy between the electroweak and 
Planck scales to exist the parameters in the Higgs sector of the standard model 
must be tuned to one part in 10"^^. New physics near the electroweak scale can 
stabilize or obviate this hierarchy. As the scale of new physics is raised above 
the TeV scale, either a new (albeit much less severe) hierarchy problem arises 
or there are problems with perturbative unitarity. Another major reason to 
consider physics beyond the standard model is in order to find a simpler or 
unified theory which offers an explanation for the standard model symmetries, 
matter content, and/or gives some framework for understanding the values of 
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the 18 standard model input parameters. 

Models of physics beyond the standard model can for the most part be 
divided into two classes, which are distinguished by their implementation of 
electroweak symmetry breaking. These are technicolor theories and supersym- 
metry. The simplest (and, it turns out, unviable) technicolor models operate 
in many respects like a scaled- up version of QCD. Such models typically en- 
counter problems related to the fermion mass hierarchy, flavor changing neutral 
currents, electroweak precision tests, and/or light pseudoscalars. To avoid this 
plethora of pitfalls, rather large, complicated, and seemingly ad hoc models 
must be invoked. 

Supersymmetry, on the other hand, does not suffer from any of the pre- 
viously mentioned problems. In fact, due to the decoupling nature of super- 
symmetric corrections, the minimal supersymmetric standard model (MSSM) 
with a heavy (~ 1 TeV) superpartner spectrum is indistinguishable from the 
standard model in all weak-scale experiments. In this light, the success of the 
standard model can be also realized as the success of the MSSM. 

The fact that the gauge coupHngs nearly unify in the MSSM El (but not 
the standard model) may be taken as indirect evidence for supersymmetric 
grand unified theories (GUTs). String theory is the only viable candidate 
for a theory of quantum gravity, and supersymmetry remains (for the most 
part) an essential element in the construction of sensible string theories. In 
some string models, low-energy supersymmetry is an unavoidable by-product. 
Gauge coupling unification can arise in string theory, with or without a GUT. 

Another nice virtue of supersymmetry is related to electroweak symmetry 
breaking. In the standard model, electroweak symmetry breaking is accom- 
plished by setting the Higgs mass parameter to the wrong-sign by hand. In 
popular supersymmetric models, however, the Higgs mass parameter has a 
positive value at an initial (high) scale. By virtue of radiative corrections the 
mass runs with the scale, and at the electroweak scale it has the wrong sign. 
This radiative breaking of electroweak symmetry is a rather generic feature of 
supersymmetric models. 

While supersymmetric theories readily explain electroweak symmetry break- 
ing, they introduce the problem of supersymmetry breaking. A priori there 
are hundreds of arbitrary new parameters in the supersymmetry-breaking La- 
grangian. There are simple models of supersymmetry breaking which greatly 
reduce the number of parameters. For example, in the simplest gravity- and 
gauge-mediated models that we will consider, there are about 5 input param- 
eters. 

With this introduction, we are motivated to study in some detail the rami- 
fications of low-energy supersymmetry. Radiative corrections play an essential 
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role in much of the diseussion of supersymmetry phenomenology. In the next 
section we will pedagogically review renormalization of the electroweak sec- 
tor of the standard model. This discussion applies to the MSSM as well. In 
the following section we discuss regularization, threshold corrections, and the 
renormalization group equations (RGEs). We apply these results when we dis- 
cuss the global fit of electroweak data in the standard model and the MSSM 
in Sec. 4. A brief summary and conclusions are provided in Sec. 5. 

2 Electroweak Renormalization 

Renormalization remains a cornerstone of particle physics. In any model of 
particle physics, whether the standard model or beyond, if radiative correc- 
tions are to be taken into account some renormalization procedure must be 
implemented. Here we review the standard technique of counterterm renor- 
malization. We will focus on the renormalization necessary to calculate the 
supersymmetric corrections to electroweak observables. For this we need not 
consider full gauge-sector renormalization (ghosts; gauge- fixing). For a more 
complete treatment see Ref. 

We first introduce the counterterms, and then discuss the renormalization 
conditions which determine them. In the following section we consider the cal- 
culation of the counterterms (and Feynman diagrams in general) by discussing 
regularization and Passarino-Veltman functions. 

The renormalization procedure is straightforward. We start in each case 
with the bare Lagrangian, which consists of bare parameters (bare masses and 
couplings) and bare fields. We write this in terms of the renormalized param- 
eters and fields. For example, the bare mass mb and coupling Cf, are replaced 
by the renormalized parameters m and e, and the associated counterterms, 6m 
and 5e, 

nib — m + 5m , Cb — e + 5e . 

The bare fields are equal to the renormalized fields multiplied by a wave- 
function renormalization factor, e.g. 

(t)b = Z^(j) . 

The zi factor is in general a matrix, mixing fields with the same quantum 
numbers into each other. After these replacements, we can separate the bare 
Lagrangian into the renormalized Lagrangian and the counterterm Lagrangian. 
The counterterm Lagrangian gives rise to new vertices, which lead to new 
graphs in each order of perturbation theory. These counterterm graphs cancel 
subdivergences in multiloop diagrams. At one-loop order, however, the coun- 
terterm vertices enter only in tree graphs. Hence, the one-loop renormalization 
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procedure is quite simple, and amounts to adding the one-loop diagram con- 
tribution to the trcc-lcvcl (bare) Lagrangian contribution. In the last step of 
the renormalization procedure we apply the renormalization conditions. The 
renormalization conditions determine the counterterms and the physical mean- 
ing of the renormalized parameters. 

For mass and wave-function renormalization, it is most natural to impose 
the on-shell renormalization conditions. In this case the renormalized mass is 
the pole mass (i.e. the experimentally measured mass). The on-shell renor- 
malization conditions consist of two parts: 

a) the renormalized mass is the real part of the pole of the propagator, (1) 

b) the real part of the residue of the pole is unity. (2) 

In the general case the inverse propagator is a matrix. The bosonic form is 

+nii(p2) -jnl^+Ui2{p^) ■ 

£,-l^p2^^| -mii-Fn2l(p2) p2_^2^^n22(p2) ••• I (3) 



Here is the external momenta-squared, rriij is the tree-level (bare) mass 
matrix and Ilij{p'^) is the one-loop self-energy. The poles are determined 
from 

Det(i?-i(p2))=0, 

and the pole masses are m| = TZe (p?). 

Scalar mass and wave-function renormalization provides the simplest ex- 
ample. The quadratic part of the bare Lagrangian is 

£ = df_,(f)ld''(l)b - ml(j)l(j)b 

= Z^d^,4>*d>'<t> - Z^{m^ + dm'Wcp . 

(4) 

So, the bare inverse propagator can be read off, and at one-loop level we simply 
add the self-energy to obtain the renormalized inverse propagator, 

D-\p-') = Z^p-" - Z^{m^ + 5m'') + n(p2) . 

The real part of the pole of the propagator is 

Tie (p2) = dm' - Z^^TZeHip^) . 

The on-shell renormalization condition requires TZe (p^) = m^. This implies 

6m^ = Z^^7een(p2) . (5) 
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At 0{a) we can set Z^^ = 1 and = in this equation, to obtain 

Sm'^ = neU{m^) . (6) 
The residue of the pole of the propagator is 

res(D(p2)) ^Z^ + n'(p2) . (7) 

The prime denotes the derivative with respect to . The second part of the 
on-shell renormahzation conditions requires that the real part of the residue 
be equal to unity, 

Z^+neU'ip^) = 1 . (8) 

Again, at 0{a) we can replace the complex pole p^ by the real part, the 
renormalized mass, so that 

= l-7een'(m2) . (9) 

Next we consider gauge-boson mass and wave-function renormahzation. 
In this case the inverse propagator has a longitudinal and transverse part. 



(10) 



Only the transverse part of the propagator contributes to physical processes. 
Correspondingly, only the transverse part of the one-loop self energy will ap- 
pear in the corrections. 

The gauge boson renormahzation is quite similar to the scalar case we just 
considered. In the standard model j-Z mixing introduces a slight complication. 
The quadratic part of the bare Lagrangian involving 7 and Z is 

K)( n n V^O , (11) 



where □ denotes the D'Alembertian. As usual we rewrite this in terms of the 
renormalized parameters and fields. 



z^, zl\(u 

Z^ZlJK U-MI-5MI 

V Z'^^U + Zl^(U -Ml- 5Ml) Zzz{n -Ml- 5Ml) 

Zf^ 





{A^'Z^') 



We have dropped some terms of 0{a^) such as (1 — Z^^)Z^2. The inverse 
propagator may be read off as 

D-^\{v^) Z,,{p^ - Ml 5Ml) + nl.ip^) ) ' ^^'> 

where 

D-zip') = D^'ip') = Z^zP^ + - - ^^^z) + n^z(p') , (13) 

and we added the transverse part of the gauge-boson self-energy, n^(p^). 

We determine the counterterms by applying the on-shell renormalization 
conditions, which in this case can be realized by|^ 

a) D-^{0)=0 (14) 
7^e^)^^(M|) = (15) 
D-'^iO) = neD-^iMl) = (16) 

b) neD-^'{Q)^TZeDz'z'{Ml) = l (17) 

These determine the counterterms 

<5M|=7^e^i^(M|) , (18) 
= l-7^e^^^(0) , (19) 
Zzz = l-nenl'ziMl) , (20) 

^'z = ~'^e-'^T^, (21) 

Z-i = ' . (22) 

Ml ^ ' 



Note that no photon mass counterterm is required since the U(l) Ward identity 
ensures n!^^(0) — 0. For the ly-boson we find 

=7^e^^^(M2,) , (23) 
Zww = 1 - 7^e^^,^.(M2,) . (24) 

The final mass and wave-function renormalization we need to consider ap- 
plies to fermion fields. Just as in the previous examples, we start by writing the 



"We use the fact that a self-energy evaluated at zero external momentum is real. 
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bare Lagrangian, then substitute for the renormahzcd fields and mass. In gen- 
eral we will have mixing among the fields, so the masses and the wave-function 
renormalization constants will be matrices. The bare kinetic Lagrangian in the 
mass basis is0 

— i+ i — 1 + 

-i^L^L {m + Sm)Z^^/jR~ iPrZ^ {m + Sm)ZLipL , 



where 4'l.r ~ (1 T 75)/2 ip- From this expression we can read off the bare 
inverse propagator, and at one-loop we simply add the self-energy contribution 
(Si, Ejj, S5) to obtain the renormalized inverse propagator|. 



•4t 



We can define Z£ ^ = 1 + 5Zl r/2. Then we have 



5Zr 



5Zl 



rriiSij 



rriidij 



SZr + 5Zl 



i>PL 
i>PR 



Pl 
Pr 



To apply the on-shell renormalization conditions it is convenient to intro- 
duce an on-shell spinor u which satisfies the equation of motion — m)u = 0. 
Recall that the on-shell renormalization conditions have two parts. The first 



''Note in the mass basis m and Sm are real and diagonal. 

■^We assumed hermiticity, which implies that S^ii and Sij^i are real and Ss/j = S, 
the fermion under consideration is unstable the effective Lagrangian is not Hermitian. 



,t 



If 
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part requires that the poles of the renormahzed propagator are at the renor- 
mahzed masses. A pole in the propagator corresponds to a vanishing column 
in the inverse propagator, when the spinor multiplies the inverse propagator 
on the right. We implement this condition by requiring 



(25) 



for all The second part of the on-shell renormalization conditions requires 
that the residues of the propagator poles are unity. Note that Eq. (|2^) only 
requires that the diagonal elements of the inverse propagator have the form 
— rrii). To ensure A = 1, we require 



1 



]/> — m 

From Eq. (pq), with i — j, we have 



(26) 



- {SZl - SZji)^,. rrii - dm^ = -miY^Ltiim^) - Sgi^ 



1 
2 



The sum and difference give 

1 



drrij — 



nii + TZe'Ssiiimf) 



5Zlu - SZjiu = Eflij(mj) - Siij(rnj) H ImY^s ii{mi ) 

rrii 



(27) 
(28) 

(29) 
(30) 



With i ^ j, we find from Eq. ( |25| ) 

^6ZLij rrij - ^jjii SZRij + S^.y (m|) nij + ^^^^(m^) = 
^SZuij rrij - ^rui SZlij + S^y (m|) + Y.sij{m'j) = 

This then determines 



(31) 
(32) 



SZ^ij 
2 



2 2 

m| - rrij 



(33) 



SZuij 
2 



(34) 



2 2 

— nij 



Eq. (pq) implies 



i) - n i i 



- [5Zl + 5Zl) J-^[SZl + SZr ] rrii 



u{mi) = , 



u{mi) = 



Taking the sum of these equations, and substituting Srrii from Eq. 
find 



1 



\ (sZl + SZl + 5Zr + SZ^J^ {i>- rui) 



where 



Note that 
1 



+ ( f{P^)i> - .f{ml)m^j + [gip') - 9{mf) 

f{p^)^^LH{p^) + ^Ru{p^) , 



u{mi) = 



i> - m 



f{p fim, )mi u{mi) = /(m, ) + 2m, / (m, ) u{mi) 



(35) 



(36) 



we 



(37) 



(38) 
(39) 
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where f'{m^) = 9/(p^)/9p^|p2^„2. Eq. (|^ imphes 



5Zr + 5Zl 



fsZL + SZl\ +2fj:L{mf)-J:R{m^)\ . (40) 
V / ii \ / a 



Substituting this in Eq. (|37|), we find 

Tie SZiti = -T,Lu{mf) - mj (^^^.(to^) + Y^'^.^^imf)) - 2m^neY.'s ^^{mj) , 
Tie 6Zru = -T,Ru{mf) - mf (^^.^(m^) + E^,^(to^)) - 2mineT,'s ^^{mf) . 



Note that the Lagrangian is invariant if we rotate both the left and right- 
handed fields 'i/'Li, '0-Ri by phases e*^'. This means that there is an arbitrari- 
ness in the diagonal wave-function renormalization constants. The on-shell 
renormalization conditions are still satisfied if we transform the SZa 's as 

SZlii e^^'SZL'ii , SZru e'-^'SZjiu ■ (41) 

No physical results depend on this arbitrariness. This freedom allows us to 
choose the SZfm to be real. We then determine the imaginary part of SZ^u 



from Eq. (30). We add it to the real part to obtain 

SZL^^ = -^Luijr^) ~ (SL,,(m2) + (42) 

2 

rrii 

The fermion renormalization constants are somewhat complicated. If one 
introduces the effective mixing matrices, which diagonalize the renormalized 
mass matrix, the wave-function renormalization takes a much simpler form. In 
particular, the off-diagonal wave-function renormalization constants (Eqs. (|3^ 
^J)) are absorbed in the construction of the effective mixing matrices. See 
Ref. d for details. 

We will be primarily interested in the e and /i lepton wave-function renor- 
malization. In the standard model, and in the supersymmetric extensions we 
will consider, XmYisa = 0. Also, there is no lepton mixing, and we can ap- 
proximate mi = 0. In this case the lepton wave- function renormalization takes 
the particularly simple form 

5Zl = -Sl(0) , SZr = -Ei^(O) . (43) 

We will use Eq. (^) repeatedly. 

It will be convenient to write the wave-function renormalization in the 
vector/axial- vector basis. Consider the wave- function renormalization of 

Ipblt^idLPL +9RPR)i^b = i^lt^iZLghPh + ZRgRPR)lp . 

It is a simple exercise to show that 



i^lt^{ZL9LPL + ZrQrPr)^ = ijj 

where 



) - {gyZa gaZv)l5 



9v = ^(5i + 9r) , 9a = ^{9L- Or) , (44) 
Z,^^{Zl + Zr), Za^\{ZL-ZR) . (45) 
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Note that = 1 + 0{a) and Za = 0{a). Hence, we introduce SZy = 1 — Z„ 
and SZa = Za- 

With these gauge-boson and fermion mass and wave-function renormal- 
izations, it is straightforward to renormalize the electric charge. The bare 
Lagrangian contains the interaction terms 



-(e -1- be) Q ipjf_,iZy - Zaj5)i> [z^. 



7Z 



(46) 
(47) 



{gvZy + gaZa) - (gvZa + gaZvhb 



^ 1 ^Iz^M 



The terms involving are (dropping terms of 0(0^)), 



Se 



1 



l + — +SZy + -6Z^^ I + —Z 



9V ryh 



QSZa + ^zl^ 



(48) 

75 }'ip ■ 



There are three form factors in the proper fermion-photon vertex, 

u2 „2 _ (49) 



-le 



ZrUp 



We can read off the tree-level vertex from Eq. (^8|) and at one-loop we simply 
add the proper vertex to obtain the renormalized vertex 



T^ik\plpj) 



7m 



(50) 



Se 



gv 



Q\l + -+SZ, + -SZ.,^ ) + ^-Zl^ + Vx(}?,p\,pI) 



QSZa + ^zl^+V2{k\plp',) 



i<^fivk ^ r /, 2 2 2\ 
^— 1^3(fc ,Pr,P2) 



The on-shell renormalization condition for the vertex requires 
r^(0,TO^,m^) = -eQ7^ 



2 ™2n /n 2 2\ 



2top 



(51) 



with the renormalized charge e equal to that measured in Thompson scattering, 
e^/47r = 1/137.065. The absence of an axial- vector contribution in Eq. ( ^l| ) is 
an automatic consequence of the axial-vector Ward identity. 



e Q SZa + eV2{Q,m^,m^) = -gaZ^^ 
11 



(52) 



There is a corresponding vector Ward identity, 

eQ SZ,+eVi{0,mlml)^ -gaZl^ . (53) 

Diagramatically, these equations correspond to (ignoring the magnetic moment 
contribution) 



+ 1/2 



+ 1/2 




= - ffa ^|^7m (1 - 75) 



In QED the right hand side vanishes {Z^^ — 0) and the usual QED Ward 
identity is rcahzcd. In SU(2)xU(l), one can consider the various contributions 
to the right hand side from the various particles in the loop. From Eq. (p2|), 

— n^2(0)/Aff . At one-loop level, all the contributions to 11^^(0) vanish 
save one, that of the W^-boson. This is true in any model with the standard 
model gauge symmetry. Hence, in the MSSM 11^^2(0) has the same value as in 
the standard model. 

If we substitute Eq. (p 
coupling as 

- eQ7^ 



into Eq. ([50|), we can write the on-shell vector 



9v - ga 



which is equal to —eQ'-fu, by the renormalization condition Eq. (b^ 



tT/ 



(0) 



gy- ga 

Qe ' 

c Ml 



(54) 
Hence, 
(55) 

(56) 



where s (c) is the sin (cos) of the weak mixing angle. It is customary to define 
Aa = 5a/a = 2(5e/e. 

Note that the vector Ward identity, Eq. (|5^) , ensures that the renormaliza- 
tion of the electric charge is universal, i.e. the flavor dependent wave- function 
and vertex corrections drop out. 

Building on these results, we can readily derive the relationships between 
the renormalized Fermi constant, the weak mixing angle, and the VF-boson 
mass. The Fermi constant is measured most accurately in muon decay. By 
considering M^-boson mediation at scales well below the VT-boson mass, the 
bare effective Lagrangian is found to have the form 



75)^r V'b7^(l - 75)C 



(57) 
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As usual this is equal to 



G,, 



V2 



^^^Zl-'^Zl^Zl^Zl^-- /^7^(1-75)V''^ V'V(1-75)V'^= . (58) 



The proper vertex function will in general contain other current bilinear prod- 
ucts besides {V - A){V - A). However, in the SM and the MSSM these terms 
are suppressed by the small electron and muon masses, so they can be safely 
ignored. The vertex function contains the following contributions: 



1^- 
e - 



.■£ii nS,-w-(o) 

'V2 M^^ 



+ 



■^2 



+ 



T r 

_i L. 



'V2 



<5box 



71 



The renormalized coupling is then 
1 + ^ + ^ {5Zi + 6Zl + 8Z2 + 5Zl') + %M 



Gii 



^vertex ~t~ ^^box 

(59) 

The muon decay rate was first calculated in the effective QED theory, where 
the photon corrections are taken into account, yielding 



1927r3 

















2tt 





(60) 



The leading two-loop correction is included by evaluating a at the scale m^. 
It is customary to take Eq. (60) as the defining equation for G^. Then, from 
(E^), we determine the counterterm (5G^, 



V2 



G^ 2 ^ ^ ^ ^ ' 



(0) 



Gp 

71 



1 + <5qed 
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^vertex ^~ ^box 

(61) 



Sqeb includes the photon correction in the effective theory, 




Hence, we have 



-5 



VB 



(62) 



where 5vb denotes the non-QED wave- function, vertex, and box-diagram cor- 
rections. 



5vB = \ {5Z1 + 5Zl + 5Zl- + 5Zl' 



+ <5box 



3ED 



(63) 



The superpartner contributions to 5vb are given in Refs. 

We now have all the ingredients necessary to determine the renormalized 
weak mixing angle. The tree-level formula applies to the bare parameters, 



Replacing the bare parameters by the renormalized ones plus the counterterms 
we find 



Expanding to first order, and using 5c? = — (5s^, we have 



2 2 

s c 



5s' 



V2G.AB 



1 + Aa 



6M^ 
-Ml 



5G, 

Ga 



(64) 



(65) 



Because of the on-shell renormalization conditions we have applied, the renor- 
malized parameters on the right hand side of this equation include the electric 
coupling measured in Thompson scattering, the physical Z-boson pole mass 
measured at LEP, and the Fermi constant determined from the measurement 
of the muon lifetime. All of these physical quantities have corresponding coun- 
terterms. Plugging in the counterterms from Eqs. (fiSh, (Eo) and (p3), we 
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have 



2 2 

s c 



1 + U^'{0) 



c 



(0) 



Ml 



- Ue 



Ml 



(66) 



'VB 



We now specify a meaning for at one-loop level. This results in a renormal- 
ization condition which determines the counterterm Ss^ . 

All the counterterms we have encountered thus far have both "finite" and 
"infinite" parts. In the calculation of Feynman diagrams one encounters di- 
vergences. A regularization method must be employed to control these di- 
vergences. Dimensional reduction is a natural regulator in supersymmetric 
theories, as will be discussed in the next section. In dimensional reduction the 
integrals are made finite by continuing spacetime to ri = 4 — 2e dimensions. 
Then, the divergences appear as poles in 1/e. We can define a renormalized 
quantity by requiring its counterterm to be purely "infinite" , i. e. proportional 
to 1/e. The renormalized quantity is referred to as a DR parameter^ and 
depends on the DR renormalization scale Q. From Eq. (|6^) we find the DR 
renormalized weak mixing angle, 



l + Ad-7^e■ 



niz(M: 



where the hat signifies i5r renormalization. We define 



Ml 



^VB 



(67) 



Afz = Ad 



Ue 



{Ml) , 11^^(0) 



Ml 



M^ 



JVB ■ 



(68) 



The leading irreducible corrections are correctly resummed by including the 
correction Afz asQ 



~2 -2 
S C 



V2G^Ml 



1 - Afz 



(69) 



In particular, the largest correction is Aa, and this form corresponds to the 
solution of the one-loop electric charge renormalization group equation. 

To calculate the M^-boson mass we introduce a DR version of the p param- 
eter, 

„ 1 M^ Ml M^ 



Mi 



Ml Ml 



''The bar in DR signifies that the 1/e pole is subtracted along with the universal integral 
artifacts ln47r — 7^;, where 7^ = 0.577216... is Euler's constant. 
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Using = M^+ Tleh^wiM^) and M| = Af| + Tie hlzi^l), we find 



Ap = Tie 



Ml 



Defining 



we can expand Eq. ((39^ to obtain 




Ana 



(71) 



(72) 



1 



(73) 



•-0 ^0 

Then, using M?y — c^pM^, we find the prediction for the M^-boson pole mass, 



clMl 



1-Ap- 



^0 ^0 



(74) 



Since the T4^-boson mass is a physical quantity, the correction is independent of 
the renormalization scheme to the order in which we are working. In particular, 
Ap + Sq/{sq — Cg)Arz is a physical observable. This combination of terms is 
gauge invariant, renormalization scale independent, and is the same whether 
hatted or unhatted. 



3 Regularization, threshold corrections, and the RGEs 

We have determined the mass, wave-function and electric charge counterterms 
by applying the on-shell renormalization conditions. These counterterms are 
comprised of two, three, and four-point diagrams corresponding to the various 
physical processes which are used to measure the input parameters. In super- 
symmetric models with the standard model gauge symmetry (e.g. the MSSM), 
the various diagrams receive extra contributions from the superpartners, but 
the forms of the counterterms remain unchanged. 

The divergences which arise in calculating the counterterms and other 
Feynman diagrams must be regulated. In order to preserve the Ward (or 
Slavnov- Taylor) identities, and to avoid spurious complications, the regulariza- 
tion method employed should preserve the symmetries of the bare Lagrangian, 
even if they are spontaneously broken. Dimensional regularization respects 
gauge symmetry, so it is widely used as a regulator of gauge theories. In di- 
mensional regularization spacetime is continued to n = 4 — 2e dimensions. For 
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71 < 4 it is clear that, for example, the integral 

J 1^ 

is ultraviolet convergent. The divergences appear as simple poles in e. 

In dimensional regularizationcl (known variously as 'naive dimensional re- 
duction', DIMR, or DREG) the dimensionality of spacetime and the (non- 
scalar) fields is continued to n-dimensions. Thus, the index of the gauge field 
Afj_ runs from to n — 1. The Dirac algebra is n-dimensional, so, for example, 
7/i7'' = 7^71.7'^ = (2 - n)ju, and so on. 

Continuing the dimensionality of the fields to n dimensions poses prob- 
lems in supersymmetric theories. Supersymmetry requires equal numbers of 
on-shell bosonic and fermionic degrees of freedom in a multiplet. Continuing 
a four dimensional gauge supermultiplet to n dimensions will spoil this equal- 
ity, so naive dimensional regularization {e.g. MS) breaks supersymmetry. One 
must add counterterms order by order in perturbation theory to restore super- 
symmetry. A new regularization method, called dimensional redujction, was 
invented for supersymmetric theories. In dimensional reduction lI (DRED), 
spacetime is continued to n dimensions while the fields remain unchanged. 
Hence, the index on the gauge field runs from to 3, and the Dirac algehm 
is 4 dimensional. There are potential problems with dimensional reductionlij. 
Ambiguities appear in calculations using dimensional reduction at high orders 
in perturbation theory. It is unknown whether these ambiguities show up in 
calculations of corrections to physical observables. Even if they did, there may 
exist a consistent prescription which gives correct and unambiguous results. 
Also, it appears thatj-the. standard mass factorization in QCD is violated in di- 
mensional reductionElltj. The situation is unclear as there has not yet been a 
complete investigation of the handling of coUinear and infrared singularities in 
massive QCD in dimensional reduction. We stress that for our purposes these 
potential problems are purely academic. It is clear that, for the superpartner 
corrections we are concerned with, dimensional reduction is a convenient and 
valid gauge symmetry and supersymmetry preserving regulator at one or two 
loop order. We will use dimensional reduction in what follows. 

We applied the on-shell renormalization conditions to determine the mass, 
wave-function, charge and Fermi constant counterterms. As mentioned in the 
previous section, it is in some cases expedient to use minimal subtraction to 
define a renormalized quantity and its counterterm. In minimal subtraction, 
naive dimensional regularization is employed, and the counterterm is defined 
to be 'purely infinite', i.e. proportional to the pole 1/e". Usually the sub- 
script 'MS' is used to signify a minimally subtracted renormalized quantity. 
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In modified minimal subtraction, the counterterm is prescribed to contain the 
pole as well as the n-dimensional integral artifacts that always appear with the 
pole, ln47r — 75. The two most common uses of modified minimal subtraction 
are in the renormalization of the weak mixing angle and the strong coupling 
constant. Such renormalized quantities are not physical observables. They are 
model dependent and renormalization scale dependent. Contributions from 
heavy particles do not decouple in such quantities (one has to implement de- 
coupling by hand). Nevertheless, they are process independent, and they can 
be useful in comparing the implications of different experimental measurements 
in a given model. 

In supersymmetry we use dimensional reduction to regulate the integrals, 
and the subtracted quantities in supersymmetric modified minimal subtraction 
are referred to as DR renormalized quantities. We use a hat to denote a DR 
renormalized quantity. 

Consider the electric charge. If we perform DR subtraction on the bare 
charge we obtain the renormalized DR charge, 

efc = e + ^ et — {Se)ao — e{Q'^) = e + Se — {de)oo (75) 

where the 00 subscript denotes the {1/ e + IhAtt — ^je) part. 



3.1 One-loop integrals 

We now discuss the calculation of Feynman diagrams. With n < 3 there is only 
one scalar function for each n-point diagram which needs to be determined to 
evaluate the one-loop Feynman integrals. All the various tensor integrals can 
be written in terms ofJiie scalar functions. The functions are as defined by 
Passarino and VeltmanE3, except we work in the metric (1, —1, —1, —1) and in 
some cases we differ by a minus sign. 
For the one-point diagram, we have 

Q'"" / (fin I 2 ^ T^M^') ■ (76) 

J (27r)" [q'' — ni'' + ie\ Idtt^ 

The renormalization scale Q is necessary in order to keep the n-dimensional 
coupling dimensionless. This integral is 

Ao{m'')^m^(l + l + \ii^) , (77) 
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where 1/e = 1/e + ln47r — ^e- The scalar two-point mtegral is defined as 

/-)4-n [ ^ _ ^ R /' 2 2 2\ 

^ J (27r)«[g'-rn2+ie][('?-p)'-"^i + *£]"16^' ^ ""^ 

(78) 

This function can be written as 

2 2 2n _ 1 , (l-a:^)m2 + xm2-a;(l-2:)p2 



Bob^m^,m^) = -- / da;ln^ '-^ — ^ ^ ^. (79) 

e 7o Q 

Such Feynman parameter integral forms are especially useful in analytically 
evaluating the Passarino-Veltman functions in special cases, such as when one 
of its arguments is zero. The explicit formula for Bq in the general case can 
be found in Ref. 

The vector two-point function is defined as 

r)4-n f _(£<! _ i p / 2 2 ™2\ 

^ J (2^)" [g2 -ml + is] [{q -pY-ml + ie] " IOtt^ ^ ' ^ ' ' 

(80) 

Bi has the Feynman parameter integral representation 

Bi{p,m(,m4)^ — - dxa;ln^ — ^ —. (81) 



2e 

The tensor two-point functions are defined as 

n^-r^fJ^ t£ (82) 

J (27r)" [g2 -ml + ie] [{q - pf -ml + ie] 

[B2iip\mlml)p^Y + B22{p\mlml)g'^-'] . 



The vector and tensor two-point functions can be written in terms of the 
scalar functions Aq and Bq as 

Bi(/,m?,m^) (83) 

= ^i^Aoiml) - Ao{ml) + {p^ + ml - ml)Bo{p^ ,ml,ml)^ , 
B2i{p\mlml) (84) 



3p 



2 



^ /'™2n 4 /_2\ I (jl I ™1 STOjA n / 2 _2 _2 



p2 



2Ao(m^)- Ao(mf)+ + ^ ^Bo(p^m^,m^) 



Ao(to^) - Ao(™^) + ( -y+ml~mi ] Bo{j}',ml,mi) 
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2 I 2 ^ 



6 1 2 



(85) 



2p2 



Aoimj) - Ao{ml) - {ml - ml)Bo{p^ ,m\,ml) 



I 2 I 2 -^2 
+ TO]^ + 7712 — — p 



The three-point scalar functions depend on the three external momenta, 
and the masses of the three internal lines. Hence, the functions take as ar- 
guments {k"^ ,p\,P2,m\,m2,m\), where mi is the mass of the internal line 
between k and pi, and m2 is the mass of the line between pi and p2- The 
scalar function is 



1 



(27r)" [q^ — m\ + ie] [{q — pi)^ — + ie] [{q ~ pi — P2Y ~ "^3 + 



167r2 



Co 



(86) 



The vector and tensor three-point function integrals are defined as 



-)4 — n 



(27r)" [q^ — m\ + ie] [{q — PiY ^ ^"2 + ^s] [{q ^ Pi ^ P2Y ^ "^3 + 



167r2 



Ciipt + Ci2p'^ 



(87) 



and 



4~ri 



(27r)" [gf2 -ml + ie] [{q - pi)^ - mj + is] [{q - Pi - P2Y - 



ie\ 



167r2 



C21 ptp'i -f C22 P^P^ + C23 + P^K} + C24 g^'' 



Analytic formulae for the C- functions can be found in Ref. [Q. The evalu- 
ation of these functions can involve large cancellations, so much care should 
be exercised in implementing algorithms. A FORTRANjjackage, FF, is available 
which accurately evaluates these integrals, and moreEj. 
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It is much more important to be familiar with the definitions and general 
properties of these functions than to be concerned with the general formula for 
them. For example, notice the general feature that the residue of the 1/e pole 
of each function is the same as the coefficient of In Q'^ . In the calculation of a 
physical quantity (such as the VF-boson mass) the 1/e divergences cancel out. 
This cancellation exactly corresponds to the cancellation of the renormalization 
scale dependence. 

It is useful to note the behavior of the Passarino-Veltman functions in the 
limit that one of the arguments is large. Looking at Bq, for example, we can 
see from Eq. ( [79| ) that in the limit M ^ m, 

1 / \ 

Bo(M2,m2,m2) = --ln-^+ const + ^ f ^ j • (89) 

This form is true for any ordering of arguments M and m. 2Bi also has the 
same form. 



Example: SUSY-QCD corrections to the top quark mass 

Given these definitions it is easy to evaluate one-loop Feynman diagrams. For 
example, consider the QCD correction to the top quark mass. The color factor 
is 

Tt,T^^ = C2{R)S,, = = ^S,, , (90) 

and the diagram is evaluated as 

\3 J (27r)» [(g2 _ m^^) + ^e] [{q - pY + le] 

^g^, [-2M(P^"^^0)+4mti3o(p^m?,0)] . (91) 

From Eq. (|||) we have 5m ^ \ [SL(m2) + Y.ii{m^)\ m + Tie Y^sim^), so the 
gluon contribution to 5mt is 

2 

<5TOt = -/%[2Bo(m?,m2,0)-Bi(m2,m2,0)]mt . (92) 
From the forms in Eqs. d?!) and (|l|), we determine 

Bo(TO^m^O) = -r +ln-^ + 2 , (93) 

Bl(m^m^O) = l + ilng + ^. (94) 
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Hence, the gluon contribution to the top quark mass in DR renormahzation is 



Srht 



3 In 



mt 



(95) 



The top quark mass is an input parameter, and as such its counterterm is 
renormahzation scale and renormahzation scheme dependent. The combination 
of counterterms and proper functions in the prediction of a physical quantity 
is renormahzation scale and renormahzation scheme independent. 

As a second example, consider the squark/gluino contribution to the top 
quark self-energy. If we ignore squark mixing, the left-handed top squark 
contributes 



83I 



d"(7 


Pr 




(27r)" 


— mi + ie 




(q — p)^ — m? + ie 



167r2 



(96) 
(97) 



so 



(98) 



Including the right-handed top squark contribution, we find the correction to 
Srut, 

(99) 



Srrit = (^Bi(mt , to?, m?^) + Bi(m^, to?, m?^)^ mt 



FromEq. (|l|), we see that in the limit to^ ^ TOg, mt, the function Bi becomes 



1 1 



B^imt,ml,mf) = - + -ln\^ +-+0 



"^?,m^ 



(100) 



After DR renormahzation the full SUSY-QCD top quark mass correction is 
then, in the limit = nni^ = ni^^ ^ TOt, m^. 



Srht 



-3 In 



In 



91 



1 / ml, mi 



mt . (101) 



The pole mass is related to the running DR mass as 

mt = mt{Q) - Srht . 



(102) 
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3.2 Threshold corrections and the RGB 



We have just calculated the SUSY-QCD correction to the top quark mass. 
The top quark mass is an input parameter in the standard model. It is not 
a prediction. Hence, the correction is Q"dependent, and in order to obtain a 
finite result in four dimensions one must prescribe a renormalization procedure. 
We used DR renormalization. 



We see that the correction (101) involves logarithms. If a logarithmic 
correction is large the validity of perturbation theory may be threatened, as 
the expansion parameter in the perturbation series becomes (a/Tr) log(M/m). 
Fortunately, these logarithmic corrections can be resummed to all orders in 
perturbation theory. The renormalization group equation serves this purpose. 

The renormalization group equation for the top quark mass can be derived 



from Eq. (101). Since the top quark pole mass is scale independent, we have 



dm. „ ^ drht dSrht 1 , "-2^^n , \ - tn\ 



dlnQ2 dlnQ2 dlnQ2 i67r2 \^3' 

(103) 

where the dots indicate gauge and Yukawa coupling contributions from the 
other (non-QCD) loops. The RGE involves the Q-dependent DR top quark 
mass and the Q-dependent DR gauge and Yukawa couplings. There are RG 
equations for these couplings as well, and together they form a set of coupled 
differential equations. There is no closed form solution to these equations, but 
it is easy to solve them numerically. 

We will now examine how the logarithmic and the non-logarithmic correc- 



tions in Eq. (101) relate to the RGE. For simplicity we will set rUg — rut- If 



TOj ^ TOt, there is no single scale where the correction ( 101 ) is small. At the 
scale rut, the squark loop contribution is proportional to \n{mi/mt)- At the 
scale Q = the gluon loop contribution involves the same large logarithm. 
The RGE resums the large logarithms. Suppose we know the top quark pole 
mass, and we want to determine the top quark running mass at scales above 



the top squark mass. To resum the logs in Eq. (101) we first apply the gluon 



correction at the scale mt- Notice there is no logarithm at this scale, so 

m{m)^{l-^-^)m. (104) 

Now we solve the RGE to evaluate this running mass at the top squark scale. 
The full RGE for the top quark mass involves all the gauge and Yukawa cou- 
plings and must be solved numerically. Here we will neglect all but the strong 
coupling so that we can solve it directly. In the effective theory below the top 
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squark mass scale, only the gluon contributes to the mt RGE. Hence, from 
Eq. (pm|), the relevant RGE is 



dmtjt) gl{t) „ 

^^ = -l^-*W' (105) 
where t = InQ^. At one-loop order, the strong coupling RGE is 

In the regime below the squark mass scale, but above the gluino mass scale, 
we have 63 = 5. The solution to this equation is 

53-^(^2) =53-'(Qi)-^ In (||) . (107) 

We can plug this into Eq. (|l05|) to solve for rht{t). We find 



rn,{Q)^rnt{m,)(^^^y' . (108) 

This expression resums the logarithmic gluonic corrections to the top quark 
mass. It is good at any scale between the top quark mass and the top squark 
mass. To determine the running top quark mass above the squark mass scale, 
we must include the top squark loop correction. To avoid a large logarithm, 
which could spoil perturbation theory, we apply the top squark correction 
at the top squark mass scale. At this scale the logarithm in the top squark 
correction vanishes. The running top quark mass just above the top squark 
mass scale is 



rhtimi) 




mt 



(109) 

From this point we can run the top quark mass to higher scales by solving 
the full MSSM RGE. Notice that above the top squark mass scale both the 
top quark mass RGE and the strong coupling RGE change, due to the squark 
contributions. 

This example illustrates how to use the RGE to resum large logarithms. 
It also demonstrates the interplay between the RGE and the logarithmic and 
non-logarithmic parts of the threshold corrections. 
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4 Z-pole precision data and calculation of the effective couplings 

We have renormalized the electroweak sector of the standard model. Since 

the MSSM has the same gauge symmetry as the standard model, this renor- 
malization also applies in the MSSM. Consider the DR weak mixing angle, for 
example. We replace the bare parameters with the renormalized ones plus a 
shift 



nab 



V2G 



to obtain 



^2 '2 
S C 



5G^){Ml 



(5M|)(1- Ad) 



(110) 



'2 '■2 
S C 



Ad 



(5M| 



1 



1- Afz 
5G,, 



G, 



n 



(0) + 2- 



(0) 



Ml 



Tie 



Ml 



(111) 
(112) 



The first 4 terms in the last expression for Afz form the oblique part of the 
correction, i.e. the gauge boson self-energy contribution. The last part, Svb, 
is the non-oblique part of the correction. In the MSSM both the oblique and 
non-oblique terms receive superpartner loop contributions. 

Applying this renormalization leads to predictions for observables in the 
MSSM which can be compared with experiment. Such comparisons can show 
which regions of supersymmetry parameter space are favored, and show that 
certain regions are inconsistent with the data. 

We already determined the prediction for the W-boson mass. We shall 
now derive expressions for other precision observables. Most of the precision 
data is due to measurements of Z-boson production and decay. We write 
the effective Z-boson-fermion-anti-fermion vertex as —ijij,{gl^ — g^^Jb) = 
—i{gj^PL + g'^ Pr). The effective couplings depend on the effective charge 
and the effective weak mixing angle as 



gf 



gf 



2s' 



2s' 



(Ts - 2Q{s-^)% gf = (Ts - Qis'^^f) , 



CI T 



9r 



(113) 



We will derive expressions for (s®*)^ and [e'^^ / s'^^d^^). First we introduce 
related observables. 
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The fermion left-right asymmetry is 

(114) 



{9f,r+{9%r (<)'+(<)' 



Note that this asymmetry is independent of the electric charge. Hence a mea- 
surement of is just a measurement of the effective weak mixing angle, 
(sj*^)^ = sin^ 9j^. In principle, each fermion has its own weak mixing angle in 
its coupling to the Z. Each fermion also couples with its own strength, which 
is measured by the partial widths, Tf. The partial widths are proportional to 

{gff)' + {gff?- 

The polarized electron beam allows SLC to measure the initial and final 
state asymmetries independently. The left-right asymmetry 

Alr = ^^1^ = Ae (115) 

is independent of the final state couplings. The polarized forward-backward 
asymmetry, on the other hand, is independent of the initial state couplings, 

,pol _ 1 {NP,F - iV-P,F) - {Np,B - N_p,b) _ 3 ^ 

- P {Np,p - N_p,p) + {Np,B - N_p,b) ~ 4 • 

P is the polarization defined by P = (Pg- ~ -Pe+)/(l ~ Pe-Pe+) and Np^p 
(Np^b) is the number of fcrmions produced in the forward (backward) direc- 
tion. SLC measures Af for f = e,i2, r, b, c. SLC also measures the heavy quark 
observables Rb^c = Tfe.c/rhad and A^^ . 

At LEP the beams are unpolarized. Hence, when the forward-backward 
asymmetry is measured, 

^..^^^^3 ^^^^^ 

a combination of the initial and final state effective weak mixing angles is 

determined. There arc other means to measure Ae independently at LEP. For 
example, a measurement of Af^ determines A^. Also, they can extract both 
Ae and At from the r-polarization measurement, 

Ar + Aef{cOSe) 2X 

'^^^''''^ = -iTam{^)' ^^"^ = TT^- ^'''^ 

Of course the various partial widths are also measured at LEP. The measure- 
ments include 

A^^^, / = e,M,r,6,c, (119) 
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-e,T 7 



(120) 



(121) 



/ 



127r reFhad 



(122) 



Chad — -T-pr 



Ml r| 




(123) 



Rb,c 




(124) 



as well as sin^ from the forward-backward charge asymmetry {Qfb)- 

The LEP and SLC collaborations do not present their results in a model- 
independent fashion. In extracting the values of observables from the data 
they assume the gauge structure of the standard model. That is, they as- 
sume the process e+e~ — + / / is mediated by standard model processes: s 
and t channel 7 and Z exchange, and box diagrams involving 7, W and Z 
exchange. The diagrams involving the photon (the s- and i-channel exchange, 
the photonic box diagrams, and the j — Z interference) are taken into account 
by the experimentalists. In other words, these contributions are subtracted 
from the data. The photonic corrections are small. On the Z-pole they are 
suppressed by a factor Tz/Mz- The genuine electroweak box diagrams are 
also suppressed, and can be neglected. In the SM and in the MSSM, the re- 
maining processes only involve the Z-fermion-antifermion vertex. Thus, by 
subtracting the photon contributions the Z-vertices are isolated. In this way 
measurements of the coupling of the Z to f f are quoted. It is good to keep in 
mind the model dependence of the quoted results. If other particles mediate 
the e+e^ ^ / / process, the quoted values of the effective Z-couphngs are in- 
correct. In the MSSM (with i?-parity), however, the only new processes which 
mediate e'^e^ —^ff are superpartner box diagrams. On the Z-pole these con- 
tributions are negligible. Hence, the effective coupling analysis is ideally suited 
to the MSSM. One needs only to compute the supersymmetric corrections to 
the effective couplings, and compare to the quoted values. 

Virtual loops involving superpartners change the effective Zff couplings. 
This is the only relevant supersymmetric correction in Z-pole data. We now 
derive the form of the Z-vertex corrections. We apply the same renormalization 
procedure as in the Sec. ||. The bare Lagrangian is of the form 



^ = - i'blfj.igvt - 9atl5) '4>bZ^ - ebQ Tphlt^i't^t ■ 



(125) 
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As usual, we apply wave-function renormalization and replace the bare cou- 
plings by the renormalized couplings plus the counterterms, 



[gv + &gv){Zv - Zaj5) - (ga + &ga){Zv - Za"t5)l5 



(126) 



(e -I- (5e)Q ■07^ (^i> - Za^i)i> 



The Z- vertex may be read off, and at one-loop order we simply add the proper 
three point function (SA(fc^,pf,p|) to obtain the renormalized vertex (we ne- 
glect the magnetic moment contribution) 



ir^(fc2,Pi,P2 
-*7m 



(127) 



gv gv ^ gv ^ gv 



-ga 



gv_ 

ga 



5ga 
ga 



1 , 5K{k\plpl) , 

- dZzz H 75 

2 ga 



At this point we can apply DR renormalization to the couplings. The 
counterterms 5gy and 5ga are then purely 'infinite' {i.e. proportional to 1/e = 
1/e -|- ln47r — 7^;) and are cancelled by the 'infinite' parts of the other terms. 
The DR couplings and cja are in tree-level relationship with e and s, which we 
have renormalized already. On the Z-pole the effective axial coupling squared 
is 



ff\2 



i2f.2 3 



1 - 



SL(0) + Sfl(0) 



gv_ 
ga 



Ei(0)~Sj^(0) 



7^e n^^(Mi)+2 



5A„(M|,0,0) 



9a 



(128) 



where we have neglected the fermion mass. From Eq. 

S2 / ;;,>2 



we have 



V2 G^Ml 1 



5Ml 
-Ml 



^/2G„M|(l + 7^e^^%M)-- 



(0) 



(129) 
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so 



l+TZe 



niz(M|) 



I]L(0) + I]fl(0) 



3« 

5a 



EL(0)-I]fl(0) 



7^e ni'^(M| 

5A„(Af|,0,0) 



(130) 



+ 2 



5a 



g'^ determines the strength of the effective couphngs. The correction to the 
effective mixing angle is determined by g'f^ /Sa^, 



9f = ^{T3~2Qsl^) 



4|Q| 



gf 
gf 



(131) 



We define k to be the ratio of s'^g to s^. The correction to k is given by 



Ak ^ 



f 1 



V4|Q|V 



- 1 



V 5q 51- 



(132) 



where Aga/ga and Agy/g^ are the corrections inside the parentheses in Eq. (127), 
in the DR scheme. Note the common Z and vector fermion wave- function renor- 
malizations cancel out, leaving 



Ak = 



1 



4|Q|V 



- 1 



El(0)-Ek(0) 



9v Ml 2gyga 
, <5A,(A/|,0,0) <5A,(Af|,0,0) 



9a 



9v 



(133) 



c II^^(M|) 
s M| 



4|Q|2s2 



- 1 



9 9 
9v - 9 a 

29v9a 



Si(0)-I]fl(0) 



Ma(A/|,0,0) <5A,(A/|,0,0) 



9a 



We obtain the effective weak mixing angle by multiplying s'^ from Eq. (69) 
with k = I + Ak. The most important supersymmetric contributions to SAy 
and SAa are listed for the Zbb coupling in Ref. |l^. Generalizing to the other 
fermion couplings is straightforward. 
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5 Comparison of supersymmetric models cmd precision data: Ax^ 
analysis 

The effective Z-boson couplings we just derived allow us to calculate the fol- 
lowing observables: 

• Z-width 




• Ratio of hadronic to leptonic width 



Thad 



• Peak hadronic cross section 



127r reFhad 



Chad 




• Ratio of heavy quark to hadronic width 



• Left-right asymmetry 



-4/ 



iaffr + iaff)' 




• Forward-backward asymmetry 



4" = l^eAf 



• Forward-backward left-right asymmetry 




-4/ 



• Effective weak mixing angle 
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The calculation of the effective Z-boson couplings allows predictions for 20 
precision observables. We can quantify the comparison of a theory with the 
measurements by forming a x^j which gives a measure of the goodness of the 
fit. We define the below. The reliability of the results of this kind of test 
increases with the number of observables. Hence, we will consider the following 
11 additional observables in the fit: 

• 3 pole masses 

The Z-mass measurement is now as precise as the measurement. We 
fix G^, but include Mz in the fit since its error is correlated with other 
measurements. The cnmbined VF-mass measurements from CDF, DO, 
UA2, and LEP 2, giveB 

Mw = 80.430 ± 0.076 GeV. 

The copibination of results of all top decay channels at CDF and DO 
gives t2l 

rat = 175 ± 5 GeV. 
We take Mz and m,t as inputs. 

• 6 low energy observables 

The weak charges of Tl and Cs hava^been measured in atomic parity 
violation experiments (APV), givingH Qw(^°^Tl) = -114.77± 1.23 ± 
3.44, and Qv^C^^^Cs) = -72.11±0.27±0.89. The deep inelastic scattering 
(DIS) experiments have produced a measipement of a linear combination 
of effective 4-Fermi operator coefficients B, k(DIS) = 0.5805 ± 0.0039. 
The neutrino scattering experiments yield determinations of leptonic 4- 
Fermi operator coefficients lil, = -M41 ± 0.015 and g""^ = -0.507 ± 
0.014. Finally, the CLEO measurement^ of B(S X^^i) yields the 90% 
confidence level interval 1 x 10"'' < B(B X.j) < 4.2 x 10"^. When 
an observable x is defined in a finite interval, there are arguments which 
suggest that a logistic transformation y — lg{x) should be performed, 
so that the new variable y is defined on the entire real axis. Including 
the transformed variable in the results in a more Gaussian shaped 
distribution in that variable. Hence, we include \g(B{B ^sj)) in 
the x^- For a variable defined on (0,1) the logistic transformation is 
lg(a;) =ln(a;/(l-x)). 

• 2 gauge couplings 

We use the constraint Aa'^^^ = 0.02817 ± 0.000620. We also include 
the constraint as{MzX-f= 0.118 ± 0.003, which we obtain by combining 
non-Z lineshape datacil. 
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These 11 observables, combined with the 20 effective couphng observables, give 
a total of 31 observables. 

5.1 Global fit in the standard model 

The measurements and experimental errors of the observables are listed in 
Table |l|, along with the best fit values of the observables in the standard 
model. Input data are as of August, 1997. The fit is performed as follows. We 
construct the by adding the square of the deviations, 

31 /^meas r.P'''''^ \ ^ 

where 0^^°^^ and 5oi is the central value and error of the measurement of the 
I'th observable, and o^"^^ is the prediction. Four of the entries in Table |^ 
are inputs, so they do not have a prediction associated with them. In the 
corresponding term in the in place of the predicted value we use the input 
value. 

We then minimize the x^ function with respect to the inputs. In the 
standard model, these are Mz, mt, Mh, Aa^^^^, and as{Mz)- In the standard 
model the x^ analysis gives rise to pertinent information about the Higgs boson 
mass. We find a standard model central value of Mh — 66 GeV, and a 95% 
confidence level upper bound of about 330 GeV. The 68% and 95% confidence 
level contours are shown in Fig. Qin the {Mh, mt) plane. 

In the MSSM, the lightest Higgs boson is constrained by the form of the 
tree-level scalar potential to be below the Z-boson mass. Radiative correc- 
tions raise_Lhe upper bound to about 125 GeV, depending on the superpartner 
spectrum c3. The lower bound on nih in supersymmetric models is 78 GeV. 
We compare the range of rn/i predicted in the MSSM with the best fit value in 
the SM in Fig. |l|. We see that the range of values predicted in the MSSM falls 
within the 68% SM Higgs boson mass range. 

5.2 Global fit in supersymmetric models 

To discuss the global fit to the data in supersymmetry we have to specify the 
soft supersymmetry breaking Lagrangian. Since we don't know how super- 
symmetry is broken, we would, in the most general analysis, consider the most 
general soft breaking terms. However, this would lead to over one hundred pa- 
rameters, with many complicated correlations needed among the parameters 
to ensure compliance with FCNC processes, rare decays, and so on. This is an 
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Figure 1: The 68 and 95% C.L. contours in the (mt, Mn) plane in the standard model. The 
range of Higgs boson mass predicted in the MSSM is also shown. 

intractable approach. Instead, we will consider specific high scale models, in 
which the soft Lagrangian takes a simple form. Besides reducing the parameter 
space to a manageable set, these models have the virtue of automatically sat- 
isfying the aforementioned constraints. We will consider two types of models, 
which are distinguished by how the supersymmetry breaking is communicated 
to the MSSM fields. 

Supersymmetric models 

In the canonical "minimal supergravity model" @, supersymmetry breaking is 
parametrized by a universal _F-term in a hidden sector. This supersymmetry 
breaking vev is communicated to the observable sector (i.e. the MSSM fields) 
by gravitational effects, giving rise to scalar masses, gaugino masses, and A- 
terms, all of order F/Mp ~ 100-1000 GeV. Since gravity is fiavor blind, it 
is assumed that these mass parameters are universal. Hence, in this minimal 
model we take as inputs a universal scalar mass, Mq, a universal gaugino mass, 
Ml/2^ and a universal trilinear scalar coupling, ^o- These inputs apply at the 
scale associated with the messenger particle, in this case the Planck scale. 
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However, for simplicity, and to avoid model dependence, we will take these 
inputs as valid at the scale where the U(l) and SU(2) gauge couplings unify. 
A/gut 10^^^ GeV. 

In simple gauge-mediated models0 supersymmetry is dynamically broken 
in a hidden sector. Interactions between the supersymmetry breaking sector 
and a standard model singlet, S, give rise to vevs in the lowest and F com- 
ponents of S. The singlet is coupled to so-called messenger fields, M and 
M, through a superpotential coupling SMM. Through this coupling the S 
vevs produce diagonal (supersymmetry conserving) and off-diagonal (super- 
symmetry breaking) entries in the (M, M) mass matrix. The messenger fields 
are charged under the standard model gauge symmetries. Hence, through the 
usual gauge interactions, the supersymmetry breaking in the messenger fields 
is communicated to the MSSM fields. This gives rise to gaugino and scalar 
masses. The masses are proportional to the gauge couplings squared times the 
ratio of the F term of S to its scalar component, A = F/S. The masses also 
depend on the representation of the messenger fields. In order to maintain the 
near unification of couplings in the MSSM, we will consider complete SU(5) 
multiplets which do not affect gauge coupling unification at one loop. We'll 
consider a messenger sector made oin^ 5+5 fields and riio 10+ 10 fields. The 
gaugino masses and the scalar masses squared are proportional to the effective 
number of messenger 5+5 fields, = + Stiiq. The trilinear A-terms are 
not generated at one-loop, so we set Aq = at the messenger scale. We will 
consider models with n^^ = 1 {i.e. = I, riio = 0, the 5+5 model) and 
nf = 3 (i.e. = 0, nw = 1, the 10 + TO model). 

In both the gauge- and gravity-mediated models the spectrum of soft 
masses is given at some high scale (either Mqut or the messenger scale M). 
In order to calculate the radiative corrections to the precision observables we 
need to determine the physical superpartner spectrum. We specify three steps 
to accomplish this. First, we need to run the soft mass parameters from the 
initial high scale down to the vicinity of the weak scale. As we describe below, 
an appropriate scale to stop the running is the squark mass scale. Second, we 
need to construct the weak scale mass matrices of the charginos, neutralinos, 
squarks, sleptons and Higgs bosons. These include the soft-breaking mass pa- 
rameters and supersymmetry conserving terms proportional to the Higgsino 
mass parameter /x, the gauge boson masses and/or the fermion masses. We 
then determine the physical superpartner mass spectrum and mixing angles 
by finding the eigenvalues and eigenvectors of the mass matrices. 
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5.3 Electroweak symmetry breaking 

In both the grajdty- and gauge-mediated models we impose electroweak sym- 
metry breakingcJ. Electroweak symmetry breaking occurs generically over the 
parameter space in both models. The RGEs of the Higgs boson soft masses 
include terms proportional to the Yukawa couplings. These terms drive the 
Higgs boson masses toward negative values. If tan (3 is not very large, the top 
quark Yukawa coupling is larger than the other Yukawa couplings, so the mass 
of the Higgs boson which couples to the top quark {H2) is driven negative. 
This is exactly what is necessary to have electroweak symmetry breaking take 
place if tan/3 > 1. In fact, because the bottom Yukawa coupling can become 
larger than the top Yukawa coupling at tan/? > 60, the mass of the Higgs 
which couples to the bottom quark (Hi) can be driven too negative at large 
tan (3 and in this case electroweak symmetry breaking does not occur. This is 
why large tan (3 values are excluded in our analysis. 

Imposing electroweak symmetry breaking allows ua to compute the heavy 
Higgs boson mass and the Higgsino mass parameter ^E. Given the soft Higgs 
masses mn^ and toj^j, we can determine the CP-odd Higgs boson DR mass 
niA, and the DR Higgsino mass parameter fj,, 

-2 - - ^'^^ . M| , (134) 



cos 2/3 

„ tan^ 3 mfr^ — m?r, 1 ~ „ ^ ^ 

^ 1 - tan2 f3 2 ^ ^ ' 

Mz is the DR Z-boson mass (A/| = A/| -I- 7?.e ^^^(A/l)), and Wlh, = rriHi - 
ti/vi. The ti are the tadpole diagram contributions. These corrections are 
necessary in order to ensure that we are at the minimum of the one-loop effec- 
tive potential. They are the effective potential corrections in the diagrammatic 
approach. The contribution to ti/vi from a particle of mass m which couples 
to Hi is of the form m? \og{ni/Q), where Q is the renormalization scale. Since 
the squarks are typically the heaviest particles, and because of the color and 
multiplicity factors, their contribution to the tadpole corrections usually domi- 
nates. The corrections to electroweak symmetry breaking can be applied most 
reliably at a scale Q where the tadpole corrections are minimized. This will 
generally be in the vicinity of the squark mass scale. This is why we stop the 
running of the soft parameters at the squark mass scale. At Q = rriq the elec- 
troweak symmetry breaking conditions can be reliably computed. At scales far 
from rriq, large logarithms make perturbation theory less trustworthy. 

We work in the convention that the gaugino masses are positive and 
can be either sign, so the sign of jjL needs to be specified. To summarize, the 
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parameter space of the minimal supergravity model is 



Mo, Ml/2, ^0, tan/?, sgn(^) , (136) 
and the gauge-mediated model parameter space is 

Af, A, tan/3, sgn(/i), nf . (137) 



5.4 Results 

We can now describe the results of the comparison of the predictions of these 
supersymmetric models with the precision observables. A more complete dis- 
cussion can be found in Ref. 129]. Choosing a point in the parameter space 



of (136) or (137), we run down the soft mass parameters to the squark mass 

scale, and check for electroweak symmetry breaking. We then determine the 

physical superpartner masses and mixings. We check if any of the superpart- 

ner masses are below the current direct search limits. If so, we disregard that 

point, and pick a new point. If all the mass bounds are satisfied, we compute 

the supersymmetric correction to each observable, and add it to the standard 

model prediction. We can then evaluate and minimize it with respect to the 

(5) 

standard model input parameters, Mz, mt, Aa^^J^^, and ag. We have to iterate 

this process since the new values of Mz, rrit, Aa\^J^^, and ag after minimization 
correspond to a different supersymmetric spectrum than the one found in the 
previous iteration. 

The supersymmetric corrections have the following characteristics. For 
very large supersymmetric masses (m ^ Mz) the supersymmetric corrections 
to the weak-scale observables decouple, and go to zero at least as fast as Mz /m. 
In the large m limit the predictions in any supersymmetric model will match 
those of the standard model, with the standard model Higgs mass equal to the 
light supersymmetric Higgs mass. As the superpartner masses become light, 
there can be relatively large corrections to the weak-scale observables, and in 
general these corrections will upset the near perfect standard model fit to the 
data shown in Table |l|. 



The oblique approximation 

Before looking at the full one-loop results, it is instructive to consider an ap- 
proximation. The full corrections can be divided up into two sets of corrections, 
those from gauge boson self-energies, and the vertex, wave-function and box 
diagram corrections. The first set, the gauge-boson self-energy corrections, 
are called the oblique corrections. They are universal in the sense that only 
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certain combinations of gauge boson self-energies appear in the corrections to 
every physical observable. In fact, in the lowest order of a derivative expansion, 
(where, for example, the derivative H' is approximated by [Il{M^) — Il{0)]/M^) 
there are only three combinations of gauge boson self-energiea, which appear. 
We consider the combinations given by Peskin and TakeuchiEHl, the S, T and 
U parameters. They are given by 



S = 
T = 

U = 



cos^ 9w {Fzz - F^^) 

"nwwio) nzz(o) 



cos 9v 



■ cos 29wF. 



sin 6w 
si-a-Ow n^z(O) 



4 sin^ e 



w 



Ml 



Ml 



Fww ^ cos^ 9\yFzz ~ sin^ 6\yFj^ 



sm26]YF^z 



4 sin^ e 



w 



a 



where Fij = {n,j{Mf ) - Il,j{0))/Mf (except F^^ = n^T,(M|)/A/|). In the 
oblique approximation the non-oblique, process specific corrections (the vertex 
corrections, wave function renormalization, and box diagrams) are neglected. 
There is a simple argument why the oblique approximation is expected to be 
a good approximation. Every superpartner (except the gluino) couples to the 
electroweak gauge bosons, so every superpartner contributes to the oblique 
corrections. Hence, the oblique corrections are enhanced by large multiplicity 
factors. For example, there are 21 sfermion contributions. In contrast, only cer- 
tain superpartners contribute to the non-oblique corrections. The non-oblique 
corrections are constrained by the specific flavor of the fermions on the external 
lines. Hence, the oblique corrections are expected in general to dominate over 
the non-oblique corrections. 

The oblique approximation gives an encapsulated view of the supersym- 
metric corrections. Rather than looking at the corrections to 30 individual 
observables, one can simply compute the corrections to the three oblique pa- 
rameters, and then consider which observables are sensitive to which oblique 
shifts. In Fig. || we show the full range of the supersymmetric corrections 
(found in a random scan of 50,000 points in parameters space) to the three 
oblique parameters in the supergravity model, versus the light chargino mass. 
We see that the corrections to T and U are always positive. The corrections 
to S and U are small, less than 0.1, while the corrections to T can be almost 
as large as 0.2. 

These oblique shifts result in shifts in the prediction of each observable. We 
can ascertain the importance of the oblique corrections in the by dividing 
the correction to each observable (due to the shift in each oblique parameter) by 
the experimental error. We show the range of the oblique corrections to various 
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Figure 2: The range of the corrections to S, T and U in the supergravity model, vs. the 
light chargino mass. 
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Figure 3: The range of the corrections to various observables due to the oblique corrections 
to 5, T and f/, in the supergravity model. The corrections are divided by the experimental 
error. 



observables in the supergravity model in Fig. ^ We see that some observables, 
such as CThad, Re and Rb receive negligible oblique corrections, while others, 
such as Tz, Ae and Mw, receive substantial oblique corrections^. 

There are two types of corrections in going from the oblique approximation 
to the full one-loop corrections. Of course the process-specific corrections must 
be added (the vertex, wave function, and box diagram corrections). One also 
needs to correct for the full gauge boson self-energies. As we have seen, some 
of the gauge boson corrections enter as derivatives of the self-energies. In the 
oblique approximation the derivatives are approximated by a difference of self- 

■^We emphasize that the size of the corrections are measured against the current experimental 
error. Fig. a does not reflect the absolute size of the corrections. 
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Figure 4: The range of the obhque (obi), non-oblique (n-obl) and full corrections to various 
observables in the supergravity model, divided by the experimental error. The non-oblique 
correction to k is neglected. 



energies. For the difference to be small, the particles in the loop should be much 
heavier than the gauge bosons. Hence, the oblique approximation should not 
work well when some superpartner masses are light, of order Mz- This is just 
what we find. In Fig. ^we compare the oblique, non-oblique and full corrections 
in the supergravity model. We show the maximum range of the corrections 
for various observables, in each case dividing by the experimental error. We 
see that for most of the observables the largest non-oblique corrections are 
about as large or larger than the largest oblique corrections. The region of 
parameter space with light superpartners has the largest corrections - and it 
is just in this region where the oblique approximation breaks down. Hence, 
when excluding regions of parameter space it is crucial to include the full one- 
loop corrections. We turn to the exclusion analysis next. We note that the 
oblique corrections to the low-energy observables k(DIS), gy', g"^, Qw{Cs) 
and Q\y(TI) are all small relative to the experimental error. Hence, the non- 
oblique corrections are bound to be small as well. We treat these observables 
in the oblique approximation in the next section. For all other observables the 
full one-loop corrections are included. 



Excluded regions of supersymmetric parameter space 

We need to specify a criteria to exclude a point in supersymmetry parameter 
space based on the test. We will use a criteria analogous to that used in 
determining the standard model Higgs boson mass limit. The 95% confidence 
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level limit on the Higgs boson mass corresponds to Ax^ = 3.84. That is, if the 
difference between the minimum value and the value at a given Higgs 
boson mass is greater than 3.84, that value of the Higgs boson mass is excluded 
by the data with at least 95% confidence. We adopt this criteria directly to the 
supersymmetric models. In each model we randomly scan over 50,000 points 
in the parameter space, and find the minimum x^, Xmin- Then, at each point 
we compute Ax^ = — Xmin- 1^ ^X^ is larger than 3.84 we consider this point 
in supersymmetry parameter space to be excluded by the data. We emphasize 
that at each point we minimize the x^ with respect to Mz, rrit, Aa^^^^ and as- 
We show the results in two dimensional subspaces. In each plot the dashed 
curve bounds the region of the parameter space which contains points which 
are excluded. The solid line bounds the region of parameter space containing 
points which are not excluded. Since we project a multidimensional parameter 
space onto two dimensions, there are regions which contain both excluded and 
unexcluded points in each figure. The regions bounded by the dashed curves 
which are outside the solid curves are absolutely excluded, independently of the 
values of other parameters. In the following we will focus on these absolutely 
excluded regions. 

In each of the following figures, we show the results for the supergravity 
model, the 5+5 gauge-mediated model, and the 10 + 10 gauge-mediated 
model, in the three panes from left to right, respectively. In each figure we 
show the results for fi > 0. As expected, the region of parameter space with 
the lightest superpartners is excluded. This is exemplified by considering the 
input parameters which directly set the scale of the superpartner spectrum. 
In Fig. ||we show the results in the {M1/2, Mq) plane (Fig. (a)), and in the 
(M, A) plane (Figs, (b) and (c)). We see that in the supergravity model we 
exclude Mq < 9 GeV and M1/2 < 105 GeV for this sign of /x. The region of 
parameter space with the smallest values of A is likewise excluded, giving the 
limit A > 36(14) TeV in the 5 + 5 (10 + 10) model. 

These limits correspond to limits on the physical superpartner spectrum. 
For example, in Fig. ^ we show the results in the (to//+ , tan /3) plane. The 
constraint from the B(i? ^ ^sl) measurement strongly excludes light charged 
Higgs masses in the ^ > case. Depending on the model, charged Higgs masses 
below 240-330 GeV are excluded. As a last example, we show the excluded 
region in the (m^, rrig) parameter space in Fig. |^|^. The squark (gluino) masses 
must be above 280 GeV (255 GeV) in the supergravity model. In the 10 -I- 10 
gauge-mediated model the squark and gluino masses are excluded below 340 
GeV. See Ref. |2^ for the excluded regions of other masses and parameters. 

•''The first two generation squark masses are nearly degenerate. Tiie squark mass plotted is 
actually mnj^ . 
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Figure 5: Excluded and non-excluded regions in the (a) supergravity model, (b) 5 + 5 
gauge-mediated model, and (c) 10 -f- 10 gauge-mediated model, in the ^o) or (M, 

A) planes, with fj, > 0. The region of parameter space where it is possible to find Ax^ > 3.84 
is bounded by the dahsed curve. The solid curve indicates the region of parameter space 
within which it is possible to find Ax^ < 3.84. The region outside the solid curve, but 
inside the dashed curve, is absolutely excluded, independently of the values of any other 
parameters. 



The results in the /z < case are typically not as strong, since the B{B ^ 
Xa'j) constraint is much weaker. Nevertheless, there are significant regions of 

parameter space excluded. In fact, for all three models under consideration we 
find the following bounds on the superpartner spectrum for either sign of /x: 

• mq> 280 GcV , 

• mg> 255 GeV , 

• m^o > 45 GeV , 

• m-+ > 195 GeV , 

• rrieL > 105 GeV , 

• m/i > 78 GeV , 

• niH > 115 GeV , 

• mA> 115 GeV , 

• mH+ > 140 GeV . 

In each of these cases the analysis gives us information about the parti- 
cle spectrum beyond that which we can currently obtain from direct particle 
searches. This provides an example of an interesting and useful application of 
the renormalization we derived in the previous section. 
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Figure 6: Excluded and non-excluded regions in the (charged Higgs mass, tan /3) plane, in the 
(a) supergravity model, (b) 5+5 gauge-mediated model, and (c) 10 + 10 gauge-mediated 
model. The curves are as described in Fig. H 
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Figure 7: Excluded and non-excluded regions in the (squark mass, gluino mass) plane, as in 

Fig. a 



6 Conclusions 

In these lectures we have reviewed the renormalization of the electroweak sector 
of the standard model. This renormalization applies directly to the MSSM. We 
reviewed the formalism which allows us to compute the supersymmetric cor- 
rections to 21 LEP and SLC observables. We discussed regularization, which 
is necessary to calculate the Feynman diagrams. In particular we discussed DR 
renormalization, which utilizes the dimensional reduction regularization appro- 
priate for supersymmetric theories. We then gave some sample calculations, 
and illustrated the relationship between the weak-scale threshold corrections 
and the renormalization group equations. 

Finally we compared the predictions for thirty-one observables in three 
supersymmetric models with the data. We saw that the regions of parameter 
space with light superpartners do not give a satisfactory global fit to the data. 
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We adopted an exclusion criteria which led to lower bounds on the various 
superpartner masses and parameters. These indirect limits are a nice comple- 
ment to the limits found in the ongoing direct production searches, and give 
us a better idea of where to, and where not to, expect supersymmetry to show 
up. 
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measurement 


SM 


pull 


Mz [GeV] 


91. 1867 ±0.0020 


91.1867 


0.0 


Tz [GeV] 


2.4948 ±0.0025 


2.4959 


-0.4 


crhad[nb] 


41.486 ±0.053 


41.478 


0.2 


Re 


20.757 ±0.056 


20.744 


0.2 




20.783 ±0.037 


20.744 


1.1 


Rf 


20.823 ±0.050 


20.789 


0.7 




0.0160 ±0.0024 


0.0163 


-0.1 




0.0163 ±0.0014 


0.0163 


0.0 




0.0192 ±0.0018 


0.0163 


1.6 


V{t) 


0.1411 ±0.0064 


0.1476 


-1.0 




0.1399 ±0.0073 


0.1476 


-1.1 


0.2322 ±0.0010 


0.2315 


0.8 


Rb 


0.2170 ±0.0009 


0.2158 


1.3 


Rc 


0.1734 ±0.0048 


0.1722 


0.2 




0.0984 ± 0.0024 


0.1035 


-2.1 




n n74i 4- n nnzis 

w.u 1 ^± zn u.uu^o 


w.u / oy 


n n 


Aim 


0.900 ±0.050 


0.935 


-0.7 


Aim 


0.650 ±0.058 


0.668 


-0.3 


Ae 


0.1548 ±0.0033 


0.1476 


2.2 


Aim) 


0.102 ±0.034 


0.148 


-1.3 


Aim) 


0.195 ±0.034 


0.148 


1.4 


Mw [GeV] 


80.430 ±0.076 


80.386 


0.6 


rat [GeV] 


175 ±5 


172 


0.6 


Qw(Cs) 


-72.11 ±0.93 


-73.11 


1.1 


Qw(Tl) 


-114.8 ±3.7 


-116.7 


0.5 


k(DIS) 


0.581 ± 0.0039 


0.583 


-0.7 




-0.041 ±0.015 


-0.0396 


-0.1 


9a 


-0.507 ±0.014 


-0.5064 


0.0 


\g{B{B ^ X«7)) 


-8.49 ± 0.45 


-7.99 


-1.1 




0.02817 ± 0.00062 


0.02817 


0.0 


as{Mz) 


0.118 ±0.003 


0.1195 


-0.5 



Table 1: Results of a global fit to the standard model. For each observable, we list the 
experimental result, the best fit result in the SM, and the pull. The pull is the difference 
between the measured value and the prediction, divided by the error. 
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